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Abstract 

Except for a few simple digital modulation techniques, derivation of average 
bit error probability over fading channels is difficult and is an involved process. In 
this letter, curve fitting technique has been employed to express bit error probabil- 
ity over AWGN of any digital modulation scheme in terms of a simple Gaussian 
function. Using this Gaussian function, a generalized closed form expression for 
computing average probability of bit error over Rayleigh fading channels has been 
derived. Excellent agreement has been found between error probabilities computed 
with our method and the rigorously calculated error probabilities of several digital 
modulation schemes. 

1 Introduction 

Bit Error Rate (BER) or bit error probability is one of the critical performance measures 
of digital communication systems. In elementary systems the channel can be modeled 
as an AWGN wherein the calculation of BER is straightforward. However, for mobile 
environments, since the signal is received via multipath with associated fading, the 
BER of AWGN channels is no longer valid. 

To compute the average BER of a given modulation scheme in fading channels, 
the corresponding BER of that modulation scheme over AWGN is averaged over the 
fading statistics of the channel. While the procedure is seemingly simple, one often 
encounters major difficulties while averaging the AWGN BER over fading channels. 

To simplify the derivation of average bit error probability over fading channels, 
there are so many approaches in the literature for approximating Q-function | Tl f2| 13] 
fe-l |5| |6| |7| and more recently |8|. In this paper, instead of direct approximation of 
Gaussian Q-function, we approximate to a high degree of accuracy, the BER itself of 
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any arbitrary digital modulation technique by a Gaussian function. This leads us to 
an exact analytical closed form expression for the average bit error probability of any 
given digital modulation technique. 



2 System model 

Assuming absence of hne of sight, the received signal r(t) in a Rayleigh fading channel 
is expressed as 

r{t) = h{t)s{t) + n{t) (1) 

where h: is the complex scaling factor corresponding to the Rayleigh multipath 
channel 

s: is the transmitted signal 

n: is the sample function of White Gaussian Noise with mean /i and variance No/2. 
Following assumptions have been made 



1 . The channel is flat 

2. The channel is randomly varying in time with each transmitted symbol getting 
multiplied by a randomly varying complex quantity h. Since h is modeling a 
Rayleigh channel, the real and imaginary parts are Gaussian distributed with 
zero mean and variance 1/2. 



3 Average bit error probability computation over the 
Rayleigh fading channel 

Let Pj^QjjMT stand for the probability of bit error for a General Digital Modulation 
Technique (GDMT) over an AWGN channel. The GDMT can be any digital modula- 
tion technique such as PSK, QPSK, QAM, FSK, ASK etc. Let P^^g'^mt denote the 
average probability of bit error for the GDMT over a Rayleigh fading channel. In the 
following section we first discuss the Gaussian function modeling of ^Vgomt which 
will be followed by derivation of an expression for J^ffg^^r based on this Gaussian 
function model . 



3.1 Approximation of bit error probability Pi^gdmt terms of 
Gaussian function 

We model the BER probability of a GDMT over an AWGN as a Gaussian function 
using simple curve fitting techniques. Towards this, we first express P^qj^j^x' §iven 
below: 

«rT = «exp[-(^^)^] (2) 

Where a, b and c are constants to be determined through simple curve fitting techniques. 
For the curve fitting, the following parameters have been used to measure the goodness 
of fit 121: 
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• Sum of Squares due to Errors (SSE) : A value closer to indicates a better fit. 

• R-square : A value closer to 1 indicates a better fit. 

• Adjusted R-square : Any value less than or equal to 1 indicates a better fit. 

• Root-Mean Square Error (RMSE) : A value closer to indicates better fit. 

Depending on the range taken by x there can be an infinite number of choices for the 
values of a, b, and c. Since values of practical interest in Eh/ No range from to only 
a few tens of dB, one needs to obtain the values of a, b, and c, only in this range of 
Eh /No and not the entire x axis. The values of constants a, b and c along with the 
values of goodness of fit for several digital modulation techniques is given in Table [T] 
andE 



Table 1: Parametric Values for Approximating Function 



Modulation Technique 


a 


b 


c 


QPSK 


0.1059 


-2.405 


4.344 


16-QAM 


0.1793 


0.3892 


8.667 


BFSK 


0.2036 


-3.056 


6.159 


BASK 


0.1059 


-2.405 


4.344 



Table 2: Curve Fitting Performance Parameter Values 



Modulation Technique 


SSE 


R-Square 


Adjusted R-Square 


RMSE 


QPSK 


2.093 X lO-*^ 


0.9998 


0.9998 


0.0002734 


16-QAM 


3.416 X lO"-* 


0.9978 


0.9978 


0.002668 


BFSK 


2.169 X 10-^ 


0.9996 


0.9996 


0.0006722 


BASK 


2.093 X IQ-*^ 


0.9998 


0.9998 


0.0002088 



It has to be specifically noted that this simple form of Gaussian function modeling 
avoids the averaging of Q-function or error function squared over the channel fading 
statistics for certain types of digital modulation schemes. 

3.2 BER Computation for GDMT over the Rayleigh fading chan- 
nel 

In (2), the value of x corresponds to Ei, /Nq or | ft, ^ £'{, /Nq for a pure AWGN or a fading 
channel respectively. Now the conditional bit error probability of a GDMT for a sample 
value of the random variable h is expressed as 

Pb\h = aexp [ - i^—^f] (3) 

with ^ — \h\'^Ei,/No. The average bit error probability of a GDMT over a fading 
channel is now computed by averaging the conditional probability given in (3) over the 
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statistical distribution of \h\'^Ei,/NQ . Note that since h has a Rayleigh distribution, 
^ will have a Chi squared distribution with degree of freedom two with pdf as given 
below LIOJ 

1 . e 



= -exp(-^) e>0 



(4) 



where 7 



7 7 

E\\h\'^]Ei,/NQ with being the expectation operator. The average bit 



error probability, 



Fading 
b,GDMT 



over the fading channel can now be computed ifTTI as 



-pFading 

^b,GDMT ~ 



(5) 



On substituting PuniO ^ndp(^) from (3) and (4) respectively in (5), we get 



ryFaa 
^b,GDMT 



ae ^ " ' —e -I di^ 
7 



On solving (6), we get 



^Fading 
b,GDMT 



27 



■ exp 



+ 



(2b _ 1\2 
V c2 -J ) 



1 + erfi- 



27' 



where erf{.) is the error function and is defined as 1121 



erf{x) — ^= / exp(— 
Jo 



(6) 



(7) 



(8) 



Hence, (7) gives a generalized closed form expression for the computation of average bit 
error probability over the Rayleigh fading channel for any digital modulation technique. It is 
worth mentioning at this stage that though in principle, the average error probability under fading 
can be computed, the process can be lengthy, and needs to be repeated when the modulation 
scheme changes, and the equations need to be re-derived. However, the novelty and utility 
of our work lies in the fact that using a simple Gaussian approximation of AWGN BER, the 
average error probability under Rayleigh fading for any given digital modulation scheme can be 
computed through a single expression, namely equation (7). 



4 Exact average bit error probability computation over 
the Rayleigh fading channel 

To verify the validity of our generalized closed form expression (7) for BER over fading chan- 
nels, we consider four common and practically used digital modulation techniques namely QPSK, 
M-QAM, M-FSK and M-ASK. The BER for these modulation techniques over AWGN can be 
found in the literature, for example in f 121. Following the approach of section III(B) and using 
Craig's fTJj expression for the Q function as given below, 

Q{x) = - / eM'—^)de (9) 
TT Jq 2 sm o 

the rigorous expressions for the average BER of the various digital modulation techniques con- 
sidered above can be arrived at and is as given below. 
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^Fading 
b,QPSK 



7+1 



where Ki = ^"l , if 2 = 



"1 = "Vir = M-i 



arctan 



01-1 



-27ra 1 



TT log2 A/ 



T^Fading 
- b.M-FSK 



M 



7 log2 M + 2 



and 



where 02 = 



ryFading 



02 



b,M-ASK — 2 

M log2 M ^ 



/327 
^27 + 2 



(10) 
(11) 



(12) 



(13) 



Exact BER for QPSK over Rayleigh Fading Channel (RFC) 
BER for QPSK over RFC using Gaussian Function 




Figure 1: Comparison of generalized and exact result for average BER computation 
over the Rayleigh fading channel for QPSK 
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Exact BER for 16-QAM over Rayleigh Fading Channel (RFC) 
BER for 16-QAM over RFC using Gaussian Function 
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Figure 2: Comparison of generalized and exact result for average BER computation 
over the Rayleigh fading channel for 16-QAM 
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Exact BER for Binary FSK over Rayleigh Fading Channel (RFC) 
BER for Binary FSK over RFC using Gaussian Function 
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Figure 3: Comparison of generalized and exact result for average BER computation 
over the Rayleigh fading channel for BFSK 
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Exact BER for Binary ASK over Rayleigh Fading Chantnel (RFC) 
BER for Binary ASK over RFC using Gaussian Function 
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Figure 4: Comparison of generalized and exact result for average BER computation 
over the Rayleigh fading channel for BASK 



5 Verification of generalized result with the exact ana- 
lytical results for various modulation techniques 

In the numerical results to be presented in this section, we have considered Ei, /No to range from 
to 50 dB. The comparison of average BER results obtained through our generalized closed 
form expression and through the exact analytical expressions for the four digital modulation 
techniques considered is as shown in Figs[T]to|4] Figure [T| and |2] shows the results for QPSK 
and 16-QAM while Fig[3]and[4]shows the results for FSK and ASK. These results clearly show 
the excellent agreement between the error probabilities computed with our method and the rig- 
orously derived exact error probabilities. 



6 Conclusion 

A closed form expression for computing the bit error probabilities of any digital modulation 
technique in a Rayleigh fading channel has been derived. Our expression for the average error 
probability does not involve any numerical methods except simple curve fitting techniques and 
is based on representing the AWGN BER in the form of a general Gaussian function. Numerical 
results show excellent agreement between error probability computed through our closed form 
expression and those obtained through rigorous derivation. 
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